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ABSTRACT 

New millimeter-wave CO and HCN observations of the host galaxies of 
infrared-excess Palomar Green quasi-stellar objects (PG QSOs) previously de- 
tected in CO are presented. These observations are designed to assess the valid- 
ity of using the infrared luminosity to estimate star formation rates of luminous 
AGN by determining the relative significance of dust-heating by young, massive 
stars and active galactic nuclei (AGN) in QSO hosts and IRAS galaxies with 
warm, AGN-like infrared colors. The analysis of these data is based, in part, on 
evidence that HCN traces high density (> 10 4 cm -3 ) molecular gas, and that the 
starburst-to-HCN luminosity ratio, L star burst /-^hcn 5 °f -f-R^-S'-detected galaxies is 
constant. The new CO data provide a confirmation of prior claims that PG QSO 
hosts have high infrared-to-CO luminosity ratios, Lm/L' co , relative to IRAS 
galaxies of comparable Ltr. Such high Li^/L' co ratios may be due to significant 
heating of dust by the QSO, or to an increased star formation efficiency in QSO 
hosts relative to the bulk of the luminous IRAS galaxy population. The HCN 
data show a similar trend, with the PG QSO host IZwl and most of the warm 
IRA S galaxies having high Ltr/Lh CN (> 1600) relative to the cool IRA S galaxy 
population for which the median (-^ir/-^hcn) CO o1 ~ 890I470. ^ the assumption is 
made that the infrared emission from cool IRA S galaxies is reprocessed light from 
embedded star-forming regions, then high values of Ltr/Lh CN are likely the result 
of dust heating by the AGN. Further, if the median ratio of L^ CN /Z/ co ~ 0.06 
observed for Seyfert galaxies and IZwl is applied to the PG QSOs not detected 
in HCN, then the derived Lir/L^ cn correspond to a stellar contribution to the 
production of L IR of ~ 7 — 39%, and star formation rates ~ 2 — 37 M yr _1 
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are derived for the QSO hosts. The corresponding values for the warm galaxies 
are ~ 10 — 100% and ~ 3 — 220 M yr -1 . Alternatively, if the far-infrared is 
adopted as the star formation component of the total infrared in cool galaxies, 
i.e., (£fir/£hcn)coo1 ~ ^starburst/^HCN^ the stellar contributions in QSO hosts 
and warm galaxies to their Lfir are up to 35% and 10% higher, respectively, 
than the percentages derived for Ltr. This raises the possibility that the Lfir 
in several of the PG QSO hosts, including IZwl, could be due entirely to dust 
heated by young, massive stars. Finally, there is no evidence that the global 
HCN emission is enhanced relative to CO in galaxies hosting luminous AGN. 

Subject headings: galaxies: active — galaxies: interacting — galaxies: ISM - 
galaxies: quasars — ISM: molecules — infrared: galaxies - 

1. Introduction 

Over the last two decades, considerable progress has been made in investigating the 
possible evolutionary connection between ultraluminous infrared galaxy mergers (ULIGs: 
i.e., Lir[8 — 1000/im] > 10 12 L ) and quasi-stellar objects (QSOs) first proposed by Sanders 
et al. (1988a,b). These studies have transformed our view of both the environments in which 
QSOs are triggered and fed, and the connection between star formation and active galactic 
nuclei (AGN) activity. 

Approximately 70% of the optically-selected Palomar-Green (PG) QSO sample have 
IRAS 1 detections consistent with thermal emission from dust (Sanders et al. 1989). In 
addition, approximately a dozen low-redshift QSOs have been detected in CO(l — > 0) to 
date (Sanders et al. 1988c; Barvainis et al. 1989; Alloin et al. 1992; Evans et al. 2001; 
Scoville et al. 2003), indicating substantial quantities of molecular gas out of which to fuel 
star formation and possibly AGN activity. 

The detection of CO(l — > 0) emission does not, in itself, imply that active star formation 
is occuring. Further, in galaxies known to harbor AGN, it is unclear what fraction of the 
infrared emission is reprocessed star light or AGN light, and thus using LpR to estimate star 
formation rates may not be justified. The present paper addresses these issues by reporting 
the results of an HCN(1 — > 0) survey of QSO previously detected in CO. HCN has a higher 
dipole moment and larger Einstein A coefficient than CO, and thus traces much denser gas 
(> 10 4 cm" 3 ). In the Galaxy, HCN emission is associated with molecular cloud cores where 
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active star formation is occuring, and it has been detected in a large number of infrared 
luminous galaxies (Gao & Solomon 2004a). Gao & Solomon (2004b) have shown that, for 
spiral and infrared luminous galaxies spanning several orders in Ljr, the infrared-to-HCN 
luminosity ratio, L^/ L' KCK , is relatively constant and has smaller deviations than Ljji/L' co 
measured over the same Li R range. Gao & Solomon conclude that HCN is a much better 
tracer of gas that ultimately gives rise to the young stellar population in these galaxies, and 
that the constancy of L IR /L' HCN is also proof that the L IR is due entirely to dust-heating 
by young, massive stars. With this as a working assumption, the data presented here will 
be used to make a comparison between the Li R / ^hcn of the QSO host sample and the 
IRAS galaxy sample. Similarities between the Li R /L' HCN ratios of the bulk of the IRAS 
galaxy population and that of QSO and IRAS galaxies with warm, AGN-like infrared colors 
are expected if the L IR of QSO hosts and warm galaxies is due to dust-heating by young, 
massive stars. 

This paper is divided into 6 sections. Section 2 is a summary of sample selection. The 
observing procedure and data reduction are presented in §3, followed by the results in §4. 
In §5, the infrared, CO, and HCN data of the PG QSO hosts are compared with that of 
both cool and warm IRAS galaxies to assess the relative contribution of star formation to 
the infrared emission of luminous AGN hosts. Section 6 is a summary of the paper. Parts 
of this discussion have been presented elsewhere (Evans 2003, 2005a, 2005b). Throughout 
this paper, we adopt an H = 75 km s _1 Mpc -1 and q = 0.5 cosmology for consistency with 
Evans et al. (2001, 2005c). 



2. Sample 

The sample was selected from Palomar-Green QSO hosts (Schmidt & Green 1983) de- 
tected in CO(1^0) to date (see Table 1). The sample is a subset of the QSO hosts im- 
aged at optical and near-infrared wavelengths by Surace, Sanders, & Evans (2001) and 
defined to have infrared-excesses (infrared-to- "big blue bump" luminosity ratio, L IR [8 — 
1000/im]/Lbbb[0.1 — 1.0/im] > 0.36) at least as high as those of ULIGs with warm, AGN- 
like infrared colors (i.e., f25fj.n1/ fm^m > 0.2: see Sanders et al. 1988b for a discussion of the 
warm ULIG sample). Within the framework of the ULIG-QSO evolutionary model, infrared- 
excess PG QSO hosts represent the stage following the warm, AGN-like ULIG phase, and 
thus their hosts should still show signatures consistent with ultraluminous, gas-rich galaxy 
mergers (tidal tails, young star clusters, molecular gas, etc.). Given this, they are also the 
most likely QSO hosts to have the highest contribution of starlight to the production of 
their infrared luminosity. The reader is referred to Surace et al. (2001) for a more detailed 
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discussion of the sample selection. 

For the purposes of the study presented here, the QSOs listed in Table 1 have three 
primary advantages. First, the CO redshifts and line widths are known, thus making it 
straight-forward to tune to the frequency of the HCN line and, in the absense of a strong 
HCN line, to measure upper limits. Second, IRAS and ISO 2 data have been published for 
these QSOs, making determinations of the infrared luminosity possible (Table 1). Third, 
because of the manner in which they are selected, they are directly comparable to the warm 
IRAS galaxy sample. 

3. Observations and Data Reduction 

Observations of the redshifted HCN emission were done at the IRAM 30m telescope 
during five observing periods between 2003 May and 2003 December (see Table 2). The 
3 and 2 millimeter receivers were used in combination with the filterbanks and an auto- 
correlator to provide bandwidths of 512 and 600 MHz, respectively. All QSO hosts with 
HCN(1 — > 0) emission redshifted to frequency accessible with the 3mm receivers were ob- 
served in HCN(1 — > 0), otherwise, HCN(2 — > 1) observations were obtained with the 2mm 
receiver. During observations, the pointing was monitored by observing standard continuum 
sources. 

Three additional sets of observations were obtained. First, all of the QSOs were observed 
in CO(l — > 0) during the 2003 December and 2006 May observing periods to verify CO 
redshifts and line-widths. Second, observations of the HCN(1 — > 0) emission from Arp 220 
and the CO(l — > 0) emission from VIIZw031 were made after every other QSO observation 
as a reality check of the 30m Telescope system performance. Finally, Arp220 was observed 
in HCN(1 — > 0) and HCN(2 — > 1) to measure the ratio of two lines (see §4). 

The data were reduced using the IRAM data reduction package CLASS. Scans were 
averaged together, and a linear baseline was subtracted from each spectrum outside of the 
velocity range of the CO and HCN line emission, or, in the cases where no HCN line was 
detected, outside the velocity range where the HCN line was expected to be. The amplitudes 
of each spectrum were converted to main beam brightness temperatures, T m b, by multiplying 
the spectrum by the ratio of the forward-to-beam efficiencies of the telescope at the specific 
frequency of the observation, then each spectrum was smoothed to ~ 20 km s _1 . Finally, 
line fluxes were measured by numerically integrating over the channels in the line profile, 
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and the line widths were measured as full width at 50% of the peak flux. The resultant 
spectra are plotted in Figures 1-3. 

4. Results 
4.1. CO 

Table 3 contains a summary of the CO emission line properties of the sample of PG 
QSO hosts. The CO luminosities are calculated via, 

L'co = 2-4 x 10 3 ( Sc ° A \ ) (£- V (1 + zr 
co VJy km s'V \MpcJ v ; 

[K km s^ 1 pc 2 ], (1) 

where Sco^v is the CO flux obtained by multiplying the line intensities in Table 3 by the 
conversion factor for an unresolved source (i.e., S/T mh = 4.95 J K+ 1 ), and D L is the lumi- 
nosity distance. The molecular gas masses, M Ha , are calculated by adopting the conversion 
factor a = M (H 2 )/L' co = 4 M (K km s" 1 pc 2 )- 1 (see Evans et al. 2001 for a detailed 
discussion of a). 

There is considerable overlap between the QSOs listed in Table 3 and recent interfero- 
metric surveys with the Owens Valley Millimeter Array (OVRO). Specifically, PG 0838+770, 
PG 1119+120, PG 1351+640, PG 1415+451, PG 1440+356, and PG 1613+658 were observed 
by Evans et al. (2001), and PG 2130+099 was observed as part of a volume-limited (z < 0.1) 
QSO CO survey by Scoville et al. (2003). In addition, previous single-dish measurements of 
CO emission from IZwl, Mrk 1014, PG 0838+770, and PG 1613+658 exist (Sanders et al. 
1988c; Barvainis et al. 1989; Alloin et al. 1992). 

The general shape and width of the CO lines are consistent with previously published 
CO spectra of these QSOs. For the QSOs with broad (> 200 km s _1 ) CO emission lines, 
the CO velocity width measured at half the maximum intensity, A-ufwhm? agree to within 
20% of prior measurements. The agreement degrades to 40% for the narrow (< 100 km s _1 ) 
CO line QSOs, the larger discrepancy resulting from the difficulty in accurately measuring 
Aw FW hm for narrow emission lines after rebinning the data. From these new measurements, 
an average CO emission line of A-ufwhm ~ 280 ± 150 km s _1 is estimated, as compared to 
the value of 300 ± 90 km s -1 measured for ULIGs (see Solomon et al. 1997 and Evans et al. 
2005c). 
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With the exception of PG 1613+658, the new measurements of main beam CO line 
fluxes, ScoA^; are on average 40% lower than the OVRO data presented in Evans et al. 
(2001). The new CO line flux measurement of PG 1613+658, which is the PG QSOs with 
the highest signal-to-noise CO detection in Evans et al. (2001), agrees to within 4% with 
the OVRO data. Finally, the CO line flux of PG 2130+099 agrees to within 10% of the 
measurement made by Scoville et al. (2003). 

While previous single-dish SqoAv measurements of IZwl and PG 1613+658 are con- 
sistent to within 20% of the measurements presented here, single-dish measurements of PG 
0838+770 and Mrk 1014 are up to 300% higher than the present measurements. The dis- 
crepancies are undoubtedly a result of calibration errors and low signal-to-noise data. 



4.2. HCN 



Table 3 also summarizes the HCN emission line properties of the QSO sample and Arp 
220. Only the PG QSO host IZwl was detected in HCN(1 -> 0); the shape and width of the 
HCN line are consistent with the CO(l — > 0) profile. For the remaining PG QSO hosts, the 
HCN line intensities were measured by numerically integrating over the channels expected to 
contain HCN emission based on the CO(l — > 0) observations. In addition, 3a upper limits 
to the HCN line intensities were derived via, 



. 3T rms Ai; FWZ i _ x 

TrcnAw < == [K km s J, (2) 

A/Aw FWZI /At; res 

where Av fwzi is the full width at zero intensity velocity of the CO emission line of the QSO 
host, and T rms is the root-mean-squared main beam temperature of the HCN spectral data 
for a velocity resolution of Af res . Line luminosities were calculated via, 



£ira , (M . 1) = 4 . 1 « xl0 .(^ r )( i ^)V.(i + ,r. 

[Kkms _1 pc 2 ]. (3) 

The HCN(l->0) lines of two of the PG QSOs (PG 0838+770 and PG 1613+658) 
are redshifted to frequencies below the range accessible with the 3mm receivers. In both 
cases, the 2mm receivers were tuned to the frequencies corresponding to the redshifted 
HCN(2^1) emission line. Both observations resulted in upper limits, which were converted 
to HCN(l^O) upper limits by adopting the £hcn(i-+o)/^hcn(2->-i) (=2-997) line ratio mea- 
sured for Arp 220 (Figure 3). 
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5. Discussion 

Approximately one quarter of the PG QSOs at z < 0.3 have been detected in CO(l— >0), 
and ~ 50 — 70% show evidence of spiral disks, merger-like tidal features, and/or luminous 
star clusters (e.g., Stockton & MacKenty 1983; Surace et al. 2001 and references therein; 
Veilleux et al. 2006; Guyon et al. 2006). The presence of star-forming molecular gas 
and luminous star clusters, combined with infrared luminosities in excess of 10 11 L , are 
compelling indications that star formation is an important component of activity in these 
galaxies, and it raises the possibility that a significant fraction of their thermal infrared 
emission may be reprocessed stellar light. 

Given such evidence, an obvious question to ask is - what are the star formation rates of 
the PG QSO hosts? As a naive first attempt, the assumption will be made that the infrared 
luminosity of the host is due entirely to dust-heating by the starburst population, and that 
r » 1 towards the star- forming regions. Doing so yields star formation rates of 

SFR ~ 1.76 x 10- 10 (^j (J^j [M yr- 1 ] ~ 20 - 420 M yr~\ (4) 

with an average star formation rate of (SFR) ~ 90 ± 125 M yr -1 and the median of 40 
M yr -1 (Table 1). Equation 4 is a modified version of Equation 4 in Kennicutt (1998), 
where the factor 5 (= 0.75 ± 0.07) is the average far infrared-to-infrared luminosity ratio, 
Lfir(40 — 500/zhi)/Lir(8 — 1000/xm), calculated for the IRA S galaxies in the Revised Bright 
Galaxy Sample (~ 600 galaxies: Sanders et al. 2003). These star formation rates are likely 
upper limits - it is difficult to imagine a geometry in which light from the AGN does not 
intersect (and thus heat) a significant fraction of dust in the host galaxy. Indeed, models 
have been developed that can account for all of infrared emission in QSO hosts by an AGN 
surrounded by a warped or clumpy torus (Sanders et al. 1989; see review in Haas et al. 
2003; Elitzur et al. 2004). 

5.1. Li R /L' co : Star Formation Efficiencies versus Dust Heating by the AGN 

To examine the importance of dust-heating by young, massive stars to the production 
of LpR in their hosts, consider first the L' co and Lir of QSO hosts relative to that of IRAS- 
detected galaxies (Figure 4: see also Evans et al. 2001, 2005c). The new L' co measurements 
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of the infrared-excess QSO hosts (Table 3) are plotted in Figure 4 3 - additional QSO data are 
compiled from Scoville et al. (2003). In order to provide a more meaningful comparison with 
the QSO hosts, the IRAS galaxies have been separated into those with warm, Seyfert-like 
(i-e., /25/im/ faofim > 0.20) and cool (/ 2 5 M m//60 M m < 0.20) infrared colors. There is significant 
evidence that most warm IRAS galaxies host AGN, but there is not clear evidence that cool 
IRAS galaxies, as a class, have luminous AGN. Given this, the cool galaxies will be treated 
as galaxies powered by starbursts. 

The new CO data show similar results to what is shown in Evans et al. (2001, 2005c) 
- the QSO hosts have low L' co for their Lt R relative to the bulk of cool IRAS galaxies. In 
other words, QSO hosts have high Z/ IR /Z/ co relative to cool IRA S galaxies with comparable 
L IR . The warm galaxies with L IR > 10 10 L show a similar trend in L IR /L co relative to 
cool IRAS sample, though not as dramatically as the QSO host population. 

To interpret the data plotted in Figure 4, consider 

-^IR -^starburst 

where L s tarburst and Lagn are the luminosities of the starburst and AGN, respectively, and e 
is the fraction of the AGN light absorbed by dust and reradiated in the thermal infrared. In 
the limit where either Lagn = or e = 0, Equation 5 is simply the star formation efficiency. 
I.e., if r >> 1 towards the starburst population which is producing stars in steady state, 
L IR /L' co is a measure of total energy output of the starburst population per unit of fuel 
available to form new stars. If e = for QSO hosts, then the high L m / L' co is an indication 
that they are producing stars extremely efficiently from the available molecular gas. Indeed, 
one would be forced to conclude that the star formation efficiency is preferentially higher in 
galaxies hosting luminous AGN. The alternative interpretation is that e > 0, and thus that 
eLAGN contributes significantly to L IR . 



3 Note that Evans et al. (2001) use Lir(1 — 1000/im) for QSO hosts instead of the commonly used 
Lir(8 — lOOOyum). The Lir(1 — lOOO^m) in Evans et al. (2001) were calculated by numerically integrating 
the spectral energy distributions in Sanders et al. (1989). Here, we use £ir(8 — lOOO^rn), which leaves out 
the significant contribution of hot dust in the 1 - 8^m range, but provides better consistency with the other 
datasets plotted. 



(5) 
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5.2. Dense Molecular Gas and Star Formation Rates 

The conclusions stated in the above discussion are only valid if CO is a robust tracer 
of star formation in these galaxies, i.e., if the CO emission is tracing most, if not all, of 
the molecular hydrogen that is actively forming stars. Given that HCN appears to be more 
closely coupled with star formation in starburst galaxies than CO (see §1), applying a similar 
analysis to HCN as done in §5.1 for CO will likely yield a more physically meaningful result. 

Figure 5 is a plot of L IR /L^ CN versus L IR of the QSO hosts, IRAS galaxies, and two 
high-redshift AGN recently detected in HCN(1 — > 0). Again, the IRA S galaxies have been 
separated into those with warm and cool infrared colors. If the dust in the cool IRA S galaxies 
are heated primarily by the starburst population (e ~ 0), then the median Li R /L' HCN of the 
cool IRAS sample will be adopted as the ratio at which L IR = L starburst . I.e., 

~ 7#7^T - 8901$ L (K km s" 1 pc 2 )" 1 , (6) 

\ L HCn/coo1 (^W«HCn) 

where o:hcn is the L' HCN -to-M^ 2 mass conversion factor and the limits on (-^ir/-^hcn) C oo1 
represent the range in which 67% (i.e., la) of the data points around the median are con- 
tained. 

The upper portion of Figure 5 (i.e., high L m /L' HCN ) is mostly populated with the warm 
galaxies and AGN hosts. Specifically, the QSO host IZwl, the two high redshift AGN 
hosts, and two of the eight warm IRAS galaxies have much higher Li R /L' HCN (> 2200) than 
the cool galaxy population (i.e., low values of ^hcn f° r their Li R ). These galaxies clearly 
have a significant contribution to their Li R from AGN activity. Further, the fact that both 
of the cool IRAS galaxies within this Li R /L^ CN range (Arp 299 and Mrk 273, for which 
Lj R /L^ CN ~ 3000) have strong X-ray signatures of embedded AGN (Delia Ceca et al. 2002; 
Xia et al. 2002; Zesas, Ward, & Murray 2003) also supports the idea that high Lir/L^cn 
are the result of significant dust-heating by AGN. Four of the eight warm galaxies have 
L m /L' UCN in the range 1600-2000, and thus overlap the high end of the bulk of the cool 
galaxy distribution. These galaxies, which include the well-studied ULIG Mrk 231, might 
have an equal mixture of AGN and starburst emission as viewed in the thermal infrared, 
however, they could also be almost entirely powered by star formation as their Li R /L^ CN is 
not outside the range of that found in cool galaxies. Finally, two of the eight warm galaxies 
have Li R /L^ CN very near the median of the cool galaxy population, illustrating the fact 
that some fraction of the warm galaxy population have "normal" Li R /L^ CN . Whether or 
not the PG QSO hosts not detected in HCN have high Li R /L^ CN similar to IZwl cannot be 
addressed with the data in Figure 5. 
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While only one of the eight PG QSO hosts was detected in HCN, and the Li R /L' HCN 
lower limits of the undetected QSO hosts are not sufficiently constraining to be conclusive, 
the Lir/L^ cn ratios of the undetected hosts can possibly be inferred by making use of 
the larger number of HCN-detected Seyferts in the IRAS galaxy sample. Figure 6 is a 
plot of L^ CN /L' co versus Lir, and is a modified version of Figure 4 in Gao & Solomon 
(2004b). Here, the IRAS galaxies have been plotted in terms of their optical emission- 
line classifications (i.e., Seyferts, LINERs 4 , and HII region-like galaxies). In this Figure, 
there is an increased spread in L' RCN /L' co values for galaxies with L IR > 10 11 L , which 
Gao & Solomon (2004b) attribute to the increase in the fraction of dense gas in luminous 
infrared galaxies corresponding to the increase in starburst activity. Eleven of the 13 (85%) 
Seyferts in Figure 6 have -^hcn/-^co < 0.1. ^ n addition, IZwl has a L' UC]<s /L' co ratio within 
the distribution of Seyfert data points at L' RCN /L' co < 0.1. Given this, it is a reasonable 
assumption that the PG QSO hosts in the present sample have low L' ncN /L ; co similar to the 
bulk of the Seyfert galaxies in Figure 6. Thus, for the PG QSO hosts, the median L' RC ^/L' co 
of Seyferts and IZwl in Figure 6, i.e. 

— — ) ~ 0.06_ 04 , (7) 

-^CO / Seyfert 

is adopted, where the uncertainties span the range of L' UCN /L' co for the 11 low Lhcn/^'co 
Seyfert galaxies. A direct estimation of Ltr/L^ cn for the QSO hosts can then be made from 
the Lm/L' co of the QSO hosts listed in Table 3. The resultant Lt R /L^ cn ratios are listed in 
Table 4, and the data are plotted in Figure 7. The stellar contributions to the production 
of Lir of PG QSOs hosts is thus 

A starburst (IR) ~ 1.76 x 10 6 ( < LlR / L HCN)eooA ( Igg/^oU 
L [K km s" 1 pc 2 ]~ 



Lir/L' ( 



CO 



7 - 39%, (8) 



and the star formation rates are 2-37 M yr -1 (median = 10 M Q yr" 1 ). Note, however, that 
given the uncertainties in Equation 7, the possibility that massive, young stars produce 60% 
of Lt R cannot be ruled out for 50% of the QSOs plotted in Figure 7. By comparison, the 
corresponding values for the warm galaxies are ~ 10 — 100% and ~ 3 — 220 M yr -1 (median 
= 20 M yr -1 ), respectively. 



4 Low Ionization Nuclear Emission Region galaxies 
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Clearly, a significant fraction of the cool IRAS galaxies do harbor AGN which may 
or may not contribute significantly to the heating of dust in their host galaxy. Thus, as a 
double-check of the results derived from Figure 7, the Figure has been replotted with the cool 
IRA S galaxies optically classified as Seyferts and LINERs 5 omitted (Figure 8). Galaxies with 
transitional spectra (e.g., galaxies classified both as LINERs and HII- region like galaxies via 
different line diagnostics) have been retained, and galaxies with no published classification 
have also been omitted. Finally, Arp 299, which has an HII region-like optical classification 
but is known via X-ray observations to harbor an AGN, has been omitted. The main 
outcome of plotting only the remaining HII region-like galaxies and transitional galaxies 
is the removal of the outliers in the L T r/L' hcn distribution. The value of (Ltr/Lh CN ) coo1 
(= 910 ± 430) remains essentially unchanged. 



5.3. Dense Molecular Gas and L FIR 

Another obvious issue to consider is whether the above results differ if the far-infrared 
emission from cool galaxies is considered as the star formation component of the total infrared 
emission, i.e., if (^fir/^hcn) CO o1 ~ ^starburst/^kcN- Figure 9 is a plot of L Fm /L' HCN versus 
Lfir of the same samples plotted in Figure 7. The use of Lfir instead of Ltr reduces the 
disparity between the PG QSO and high-z AGN hosts data points relative to the IRAS 
galaxy data points, and shows the possibility that the far-infrared emission in several PG 
QSO hosts could be due almost entirely to dust heated by young, massive stars (see Table 
4). In particular, IZwl has a Lfir/^hcn ver y near the median value of the cool IRAS galaxy 
population. In the case of the warm galaxies, the median (-^FiR/-^HCN) W arm (~ 1220+sio) 
is 33% lower than (LiR/L' HCN ) warm , however, the median value of the cool IRAS galaxies, 
(L FIR /Z^ CN ) cool ~ 660+^, is 25% lower than (L T r/L' hcn ) coo1 . Thus, the relative spacing 
of the warm and cool galaxy data points are only marginally changed by using Lfir. As a 
result, the contribution of dust-heating by young, massive stars to the production of Lfir 
for the PG QSO hosts and warm galaxies increases by 1-35% and 1-12%, respectively, over 
the contribution derived for Ltr (see Table 4). These results reflect the fact that the average 
Lfir/Lfr for the cool, warm, and PG QSO galaxies in Figure 8 are 0.77 ± 0.05, 0.64 ± 0.11, 
and 0.44 ± 0.10, respectively. 



5 Note that, while some of the LINERs are likely primarily starburst galaxies, some of them have strong 
evidence of embedded AGN. Thus, the conservative approach of omitting all of the LINERs has been taken 
for this exercise. 
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5.4. Comparison with Recent Estimates of Star Formation Rates of PG QSO 

Hosts 

The general conclusion that the PG QSO host star formation rates in Equation 4 may, 
in some cases, be overestimates is consistent with a recent study of the star formation rates 
of QSO hosts by Ho (2005). In this study, star formation rates were derived by first assuming 
that the [O II] A3727 / [O III] A5007 line ratio in QSOs is the same as the constant value 
measured for the narrow-line regions of Seyfert galaxies, then attributing excesses in the 
ratio to star formation within the host. Doing so yields [O II] A3727-derived star formation 
rates of ~ 0.6 — 20 M Q yr -1 - i.e., significantly lower than what would be derived based 
solely on the assumption that the Li^ ~ -^starburst > and is consistent with the low end of the 
range in star formation rates estimated for individual hosts in Table 4. Note, however, that 
the [O II] A3727 emission line is subject to significant extinction by dust, and thus the Ho 
(2005) estimates are likely lower limits. 

Another study (Schweitzer et al. 2006) addresses the issue of the nature of far-infrared 
emission in PG QSO hosts by applying the 7.7 /zm PAH-to-far infrared luminosity ratio 
(Lpah/^fir) versus [Ne II] 12.5/iin-to-far infrared luminosity ratio (L[ NeII ]/^Fm) diagnostic 
of starburst galaxies to the QSO host sample. Starburst galaxies are observed to have strong 
PAH and weak [Ne II] relative to AGN-dominated hosts because the AGN efficiently destroys 
PAH and is able to produce extended semi-ionized regions in which neon is collisionally 
excited. Schweitzer et al. conclude that 30% or more of the far-infrared emission to due to 
dust heating by young, massive stars, relative to the median percentage of 427^5% estimated 
from the present analysis of the HCN data. 

5.5. The Effects of the AGN on HCN Emission 

Throughout this paper, the assumption has been made that HCN is collisionally excited 
by molecular hydrogen. Recent papers have challenged this claim by stating that, in galaxies 
harboring AGN, X-ray emission from the AGN enhances HCN emission (e.g., Kohno et al. 
2003; Gracia-Carpio et al. 2006). The high Lir/L^ cn , and thus the low L' KC]SS , °f mos t of 
the PG QSO hosts, high- z AGN hosts, and warm galaxies relative to cool IRAS galaxies 
(Figure 7) clearly contradicts this hypothesis. 

To examine this issue another way, consider again the plot of L' nc ^/L' co versus Ljr 
of the QSO hosts and IRAS galaxies (Figure 6). In this Figure, IZwl and most of the 
Seyfert galaxies have -^hcn/-^co < 0.1, consistent with the distribution of the HII region- 
like galaxies. Two of the Seyfert galaxies have L' nc¥! /L' co > 0.1, however, they are still 
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within the distribution of ^hcn/^co < 0-1 °f HII region-like and LINER galaxies. Thus, 
there is no evidence that the global HCN emission is enhanced relative to CO by the presence 
of a powerful AGN. 

6. Conclusions 

New CO and HCN observations of a sample of z < 0.17 optically-selected PG QSO hosts 
with infrared excesses are presented. These data have been compared with similar data of 
IRAS-detected galaxies to assess the validity of using the infrared luminosity to estimate 
star formation rates in luminous AGN hosts. The following conclusions are reached: 

1. The new CO measurements confirm the results presented in Evans et al. (2001) - 
that QSO host as a class may have low values of L' co for their Li R relative to IRAS galaxies, 
or relatively high values of Lir/L' co . Such high Lir/L' co ratios may be evidence that AGN 
contribute significantly to the heating of dust in their host galaxies, or that they have high 
star formation efficiencies. 

2. A comparison of the QSO hosts and IRAS galaxies surveyed in HCN(l^O) to date 
show the upper distribution of L IR /L' HCN (> 1600) to be mostly populated by warm IRAS 
galaxies and AGN hosts. This distribution overlaps at the low end with cool IRAS galaxies 
for which (Lm/ L'ucn) C oo1 ~ 890^479 L (K km s _1 pc 2 ) -1 . There are also 2 (out of 8) warm 
galaxies with Lir/L' hcn comparable to the median Lir/L^ cn of cool galaxies. An enhanced 
Lir/L' hcn relative to cool IRA S galaxies thus appears to be an indication that an AGN is 
contributing significantly to heating the dust. However, the possibility that the dust in most 
of the QSO hosts and warm galaxies are heated primarily by young, massive stars cannot be 
entirely ruled out. If the median ratio of L' RCN /L' co ~ 0.06 observed for Seyfert galaxies and 
the PG QSO host IZwl is applied to the PG QSOs not detected in HCN, then the derived 
Z/ir/ L' hcn correspond to a stellar contribution to the production of Lir of ~ 7 — 35% and 
star formation rates ~ 2 — 37 M Q yr -1 are derived for the QSO hosts. The corresponding 
values for the warm galaxies are ~ 10 — 100% and ~ 3 — 220 M yr _1 . 

3. The average (Lfir/^ir) of the cool, warm, and PG QSO galaxies considered in this 
study are 0.77±0.05, 0.64±0.11, and 0.44±0.10, respectively. Thus, the disparity between 
the QSO (and high-redshift AGN) hosts and the cool galaxies is diminished by comparing 
the Lpin of QSO hosts and cool galaxies. If the far-infrared is adopted as the star formation 
component of the total infrared in cool galaxies, i.e., (-^fir/-^hcn) C oo1 ~ ^starburst/^HCN) then 
five of the nine QSO hosts, including IZwl, potentially have all of their L F ir generated by 
dust heated by massive stars. The stellar contributions to L FIR in QSO hosts are up to 35% 
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higher than those derived for Lir, but the warm galaxies have contributions only 10% higher 
than those derived for L IR . 

4. The PG QSO host IZwl and the IRAS galaxies with Seyfert optical classifications 
have L' HCN /L' co distributions comparable to Hll-region like galaxies. It is thus unlikely that 
the global HCN emission in galaxies harboring luminous AGN is enhanced relative to CO 
by radiation from the AGN. 

A HCN survey of a more statistically complete sample of IRAS galaxies and PG QSO 
hosts is required to assess the robustness of many of these results. The apparent faintness 
of HCN emission in PG QSO hosts, and in many warm galaxies, necessitates the use of 
a larger millimeter-wave telescope. Such observations may be possible with the soon-to- 
be commissioned 50m-diameter Large Millimeter Telescope (LMT), or the lOOm-diameter 
Green Bank Telescope (GBT). Frequencies as high as 90 GHz should be accessible with the 
GBT in a few years. 

We thank the telescope operators and staff of the IRAM 30m telescope for their support 
both during and after the observations were obtained, and the anonymous referee for critical 
comments that greatly improved the accuracy and clarity of the paper. ASE also thanks J. 
Mazzarella and J. Mulchaey for useful discussions and assistance. ASE was supported by 
NSF grant AST 02-06262. This research has made use of the NASA/IPAC Extragalactic 
Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronautics and Space Administration. 
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Figure Captions 

Fig. 1. — IRAM 30m spectra of CO(l — > 0) emission from nine PG QSO hosts with infrared- 
excesses. The intensity scale is in units of main beam brightness temperature. A linear 
baseline has been subtracted from each spectrum; the baseline subtraction is performed 
outside of the velocity range of the emission lines. The zero velocity corresponds to the 
redshifts listed in Table 1. 

Fig. 2.— IRAM 30m spectra of HCN(1 -> 0) and HCN(2 -> 1) emission from the QSO 
host sample. The intensity scale is in units of main beam brightness temperature. A linear 
baseline has been subtracted from each spectrum; the baseline subtraction is performed 
outside of the velocity range where the emission line is expected to be based on the position 
and velocity width of CO(l — > 0) emission from these QSO hosts. 

Fig. 3.— IRAM 30m spectra of HCN(1 -> 0) and HCN(2 -> 1) emission from the ultra- 
luminous infrared galaxy Arp 220. The intensity scale is in units of main beam brightness 
temperature. A linear baseline has been subtracted from each spectrum; the baseline sub- 
traction is performed outside of the velocity range of the emission lines. The zero velocity 
corresponds to a redshift of 0.0181. 

Fig. 4. — L' co vs. Lir of the PG QSO hosts and a sample of i7L4,5-detected galaxies. 
References for the data are as follows: infrared galaxies, Mazzarella et al. (1993a) and 
Solomon et al. (1997); PG QSO hosts, this paper and Scoville et al. (2003). 

Fig. 5. — L IR /L^ CN vs. L m of the PG QSO hosts, a sample of local IRAS galaxies, and 
two high-z AGN (FSC 10214+4724 and the Cloverleaf Quasar). Filled squares with arrows 
represent 3a lower Ltr/L' hcn limits of the PG QSOs not detected in HCN. The dashed line 
is L IR /L' HCN = 890 L Q (K km s _1 pc 2 ) -1 (see §5.1) - the error bar associated with the 
dashed line encloses 67% (la) of the cool galaxy data points. With the exception of NGC 
4418 (Imanishi et al. 2004), Mrk 273 and IR 17208-0014 (Imanishi, Nakanishi, & Kohno 
2006), and Arp 220 (this paper), the IRAS galaxy HCN data are compiled from Gao & 
Solomon (2004a,b). The F10214+4724 and Cloverleaf data are from Vanden Bout, Solomon, 
& Maddalena (2004) and Solomon et al. (2003), respectively. 

Fig. 6. — The -^hcn/-^co versus -^m of IZwl and IRAS galaxies. The IRA S galaxy data are 
plotted by their optical emission-line classification. The warm IRAS galaxies are encircled. 



This preprint was prepared with the AAS IATgX macros v5.2. 
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Fig. 7. — Same as Figure 5, except the estimated Lir/ L' ncN (see §5.2) of the PG QSO hosts 
not detected in HCN are plotted and the y-axis upper limit has been increased. The vertical 
bars associated with the PG QSO host data points cover the range in L IR /L^ CN that result 
from the uncertainties in Equation 7. 

Fig. 8. — Same as Figure 7, except only the cool IRAS galaxies optically classified as HII- 
region like galaxies or transition galaxies are plotted (Arp 299 has also been omitted - 
see text). The warm galaxies and PG QSO host data points remain unchanged. Optical 
emission-line classifications are obtained from Veilleux et al. (1995), Ho, Filippenko, & Sar- 
gent (1977), Mazzarella & Boroson (1993b), and the NASA Extragalactic database (NED). 

Fig. 9. — Lfir/^hcn vs - -^ fi R °f the QSO hosts, a sample of local IRAS galaxies, 
and two high-z AGN (FSC 10214+4724 and the Cloverleaf Quasar). The dashed line is 
Lfir/£' hcn = 660 L (K km s" 1 pc 2 ) -1 (see §5.3) - the error bar associated with the dashed 
line encloses 67% (la) of the cool galaxy data points. 
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Table 1. Source List 





Coordinates (J2000.0) 


D L a 




T a 
-^IR 


T a 
^FIR 


SFR b 


Source 


R.A. 


Decl. 


(Mpc) 


zco 


(xlO u L ) 


(xlO n L ) 


n r 1 

M yr- 1 


pp nnc;n_i_i9/i T7-nri 

trKji UUDU-rlZ4 — IZjWI 


00:53:34.94 


+ 12:41:36.20 




U.UOl 




Z.o 


yo 


PG 0157+001 = Mrk 1014 


01:59:50.21 


+00:23:40.62 


680 


0.163 


32 


21 


420 


PG 0838+770 


08:44:45.36 


+76:53:09.20 


540 


0.131 


2.9 


1.4 


40 


PG 1119+120 = Mrk 734 


11:21:47.15 


+ 11:44:19.00 


200 


0.050 


1.1 


0.44 


15 


PG 1351+640 


13:53:15.83 


+63:45:45.60 


360 


0.088 


5.6 


2.1 


75 


PG 1415+451 


14:17:00.84 


+44:56:06.50 


470 


0.114 


1.6 


0.66 


20 


PG 1440+356 = Mrk 478 


14:42:07.48 


+35:26:23.10 


320 


0.078 


2.8 


1.4 


40 


PG 1613+658 = Mrk 876 


16:13:57.21 


+65:43:10.60 


530 


0.129 


7.6 


3.9 


100 


PG 2130+099 = UGC 11763 


21:32:27.77 


+ 10:08:19.50 


250 


0.063 


2.1 


0.62 


30 


Average 














90 ± 125 


Median 














40 



a Calculated using the prescription in Sanders & Mirabel (1996), and assuming H = 75 km s^ 1 Mpc -1 and 
g = 0.5. The Li R of PG 1119+120 and PG 2130+099 are calculated using the flux densities published by Haas et 
al. (2003). For the remaining QSO host galaxies, the flux densities published by Sanders et al. (1989) are used. 

b Calculated by assuming that L m = L sta rburst (see §5). 



Table 2. Journal of Observations 



Date Time T sys 



Source Transition (mmm-yy) (hr) (K) 



PG 0050+124 CO(l-> 0) 

HCN(1 -> 

HCN(1 -> 
PG 0157+001 CO(l-> 0) 

HCN(1 -> 
PG 0838+770 CO(l-> 0) 

HCN(2 -> 1 
PG 1119+120 CO(1->0) 

HCN(1 -> 
PG 1351+640 CO(l-> 0) 

HCN(1 -»• 
PG 1415+451 CO(l-> 0) 

HCN(1 -> 
PG 1440+356 CO(l-> 0) 

HCN(1 -> 
PG 1613+658 CO(l-> 0) 

HCN(2 -> 1 
PG 2130+099 CO(l-> 0) 



Dec-03 


40 


161 


May-03 


527 


102 


Sep-03 


215 


111 


Dec-03 


50 


111 


Dec-03 


70 


260 


Dec-03 


45 


111 


May-03 


220 


217 


Dec-03 


67 


95 


Jun-03 


301 


108 


May-06 


98 


146 


Jul-03 


300 


111 


Dec-03 


50 


106 


Dec-03 


205 


116 


Dec-03 


80 


126 


Dec-03 


210 


120 


Dec-03 


35 


126 


Dec-03 


125 


235 


Dec-03 


140 


157 



Table 3. Emission Line Properties 



Source 


zco 


line 


A?JFWHM 
(km s" 1 ) 


T mb Ai> 
(K km s" 1 ) 


( Jv km s — M 


jj a,b 
line 

(K km s- 2 pc 2 ) 


M(H 2 ) C 


L d 


j d,e 
^IR 

L l 

HCN(1 — >0) 


PG 0050+124 


0.061 


CO(1->0) 


370 


6.0+0.2 


30+1.1 


4.2 X 


10 9 


1.7 x 10 10 


170 






0.061 


HCN(l^O) 


290 


0.27+0.03 


1.3+0.1 


3.2 X 


10 8 




2260 




PG 0157+001 


0.163 


CO(1->0) 


270 


1.1+0.1 


5.5+0.5 


5.2 x 


10 9 


2.1 x 10 10 


610 








HCN(l^O) 




-0.12 ± 0.16 


-0.60 + 0.80 


< 9.6 x 


10 s 




>3300 


>1100 


PG 0838+770 


0.132 


CO(1->0) 


60 


0.50+0.07 


2.5+0.4 


1.5 x 


10 9 


6.1 x 10 9 


190 








TJ (~\ "NT / O . -1 \ 

HOIN [2 — >L ) 




n 1 o _L n c\a 
U.lz ± U.Ub 


U.bU ± U.ol 


< 1.8 x 


10 s 




>1590 


>530 


PG 1119+120 


0.050 


CO(l-vO) 


220 


0.54±0.09 


2.7±0.4 


2 5 x 


-LU 


1.0 x 10 9 


450 








HCN(l^O) 




-0.01 ±0.04 


-0.05 ±0.19 


< 9.2 x 


10 7 




>1230 




r{j 1O01+D4U 


U.Uoo 


— >U ) 




n c^/i-Ln no 
U.04±U.Uy 


z. / ±U.O 


7.7 x 


10 s 


o.l X 1U 


1 oU 








HCN(l^O) 




0.09 ± 0.08 


0.46 ± 0.40 


< 5.5 x 


10 s 




>1000 




PG 1415+451 


0.114 


CO(l— >0) 


90 


0.42+0.06 


2.1+0.3 


9.8 x 


10 s 


3.9 x 10 9 


160 








HCN(l^O) 




-0.02 ± 0.04 


-0.12 + 0.20 


< 3.8 x 


10 s 




>420 




PG 1440+356 


0.078 


CO(l-vO) 


310 


1.3+0.1 


6.6+0.6 


1.5 x 


10 9 


2.0 x 10 9 


190 








HCN(l^O) 




0.09 ± 0.08 


0.43 + 0.40 


< 3.9 x 


10 s 




>720 




PG 1613+658 


0.129 


CO(1->0) 


400 


1.6+0.1 


8.0+0.6 


4.8 x 


10 9 


1.9 x 10 10 


160 








HCN(2^1) 




-0.17 + 0.13 


-0.83 + 0.63 


< 4.1 x 


10 s 




>1850 


>620 


PG 2130+099 


0.063 


CO(l-vO) 


530 


0.78+0.1 


3.9+0.5 


5.6 x 


10 s 


2.2 x 10 9 


380 




Arp 220 


0.0181 


HCN(l^O) 


550 


9.0 + 0.8 


45 + 3 


9.6 x 


10 s 




1450 








HCN(2->1) 


550 


12 + 1.4 


60 + 7 


3.2 x 


10 s 








PG QSO average (CO) 






260+150 

















a 3cr upper limit. 

b Calculated assuming Ho = 75 km s _1 Mpc -1 and qo = 0.5. 
c Calculated assuming a = 4Mq (K km s _1 pc 2 ) -1 . 
d In units of L Q (K km s" 1 pc 2 )" 1 . 

ei HCN(2^i)/^HCN(i^0) conversion is achieved by adopting the ^hcn(1->o)/^hcn(2- > i) (= 2 - 997 ) °f Arp 220. Sec §4.2. 
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Table 4. Infrared-to-HCN Luminosity Ratio Comparisons and Adjusted Star Formation 

Rates 



Source 


Type 


Class a 


-km b 


^FIR b 


It it ' \ 

( L IR/ L HCN/ coo | 


It it' \ 


SFR d 


L nc.N 


r 1 
HON 


Lm/L'nr,N 




median 


cool 




890±*« 


660+ 370 

DDU -330 








median 


warm 




1850™ 


122011?" 


0.48 ± 0.25 


0.54 + 0.28 




median 


QSO 




qi r n +6300 


1 575+ 4720 
lo,o -530 




n 4O+0.56 
u -+^-0.35 




PG 1351+640 


QSO 


Seyfert 1 


1 91 %+™« 


4540+ yuyu 
404u -1520 


u - u '-0.06 


15+ U - Iy c 
u - io -0.12 


5 +Y 

°-5 


PG 0157+001 


QSO 


Seyfert 1 


10120™ 


6640+ 13280 


09+ 011c 
u - uy -0.07 


io+°- 13c 




PG 1119+120 


QSO 


Seyfert 1 


7 ro n +15070 




12+ 015c 


n 99 +0. 29c 


2± 2 


Cloverleaf Quasar 


high-z AGN 


Seyfert 1 


7700 


1700 


0.12 ± 0.06 


0.39 + 0.20 


375 ± 190 


PG 2130+099 


QSO 


Seyfert 1 


6250+ 12500 


1840+ 3 ^ 


n 14+o.iSc 


n q^+0.48c 
U - OD -0.30 


4+5 
4 -3 


IRAS FSC10214+4724 


high-z AGN 


Seyfert 2 


4800 


2600 


0.19 ± 0.10 


0.25 + 0.13 


223 ± 110 


NGC 4418 


warm 


LINER 


4610 


3830 


0.19 ± 0.10 


0.17 + 0.09 


3 + 2 


PG 1440+356 


QSO 


Seyfert 1 


,i c n +6300 


1580±Hg° 


0.28t°' 3 5c 


n 49+0-56c 




PG 0838+770 


QSO 


Seyfert 1 


oidu_ 1050 


1520t 3 5 ° u 


0.28i8'g c 


n 4 o+0.58c 


n± 9 3 


PG 1415+451 


QSO 


Seyfert 1 


2720±^ 30 


1120t 22 ° 


0-33±°;« c 


U-59_o.4 9 




PG 1613+658 


QSO 


Seyfert 1 


2630™ 


1380 2 _ 7 4 6 6 ° 


n 4 +0.42c 


0.48_ 40 


o 4 +42 


NGC 1614 


warm 


HII/LINER 


3090 


2090 


0.29 + 0.15 


0.32 + 0.17 


15 + 7 


PG 0050+124 = IZwl 


QSO 


Seyfert 1 


2260 


880 


0.39 ± 0.20 


0.75 + 0.39 


37+ 19 


IRAS FSC05189-2524 


Seyfcrt 2 


warm 


1900 


1260 


0.47 ± 0.24 


0.52 + 0.28 


73 ± 37 


NGC 7469 


warm 


Seyfert 1 


1860 


1190 


0.48 ± 0.25 


0.55 + 0.29 


26 ± 13 


NGC 3034 = M82 


warm 


HII 


1850 


1270 


0.48 ± 0.25 


0.52 + 0.27 


3 + 2 


Mrk 231 


warm 


Seyfert 1 


1630 


1030 


0.55 ± 0.28 


0.64 + 0.34 


219 ± 112 


NGC 1068 


warm 


Seyfert 2 


780 


320 


1.1 + 0.59 


2.0+ 1.1 


42 ± 21 


NGC 7479 


warm 


Seyfert 2 


660 


410 


1.4 ± 0.70 


1.6 ± 0.83 


13 + 7 



a References for optical spectral classification: Veillcux ct al. (1995) and Ho, Filippenko, & Sargent (1997). 

b In units of L (K km s" 1 pc 2 )" 1 . 

c Calculated by adopting L' HCN /L' co = 0.06 (see §5.1). 

d Star formation rate in Mq yr _1 . 



